Systematic first-principles study of impurity hybridization in NiAl 
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We have performed a systematic first-principles computational study of the effects of impurity 
atoms (boron, carbon, nitrogen, oxygen, silicon, phosporus, and sulfur) on the orbital hybridization 
and bonding properties in the intermetallic alloy NiAl using a full-potential linear muffin-tin orbital 
method. The matrix elements in momentum space were used to calculate real-space properties: 
onsite parameters, partial densities of states, and local charges. In impurity atoms that are empiri- 
cally known to be embrittler (N and O) we found that the 2s orbital is bound to the impurity and 
therefore does not participate in the covalent bonding. In contrast, the corresponding 2s orbital 
is found to be delocalized in the cohesion enhancers (B and C). Each of these impurity atoms is 
found to acquire a net negative local charge in NiAl irrespective of whether they sit in the Ni or Al 
site. The embrittler therefore reduces the total number of electrons available for covalent bonding 
by removing some of the electrons from the neighboring Ni or Al atoms and localizing them at the 
impurity site. We show that these correlations also hold for silicon, phosporus, and sulfur. 
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I. INTRODUCTION 

The development of better nickel-based superalloys has 
paced the construction of larger, more powerful, and 
more fuel efficient aircraft and industrial gas turbines. 
jl],^| Aluminum is the most important alloying element in 
nickel, for both strength and oxidation resistance. jl],|]]U 
Pure NiAl, which crystallizes in the B2 structure, has 
low density, high melting temperature of 1638°C (melt- 
ing temperature for the fee nickel is 1455°C), and good 
electrical and thermal conductivity. Its practical ap- 
plication, however, is limited by poor toughness and dam- 
age tolerance at room temperature j^] and brittle grain- 
boundary fracture at ambient and elevated temperature. 
H The strength and other properties of NiAl can be mod- 
ified by adding various impurity atoms. Typical mod- 
ern nickel-base superalloys contain eight or more differ- 
ent elements, each with specific functions with respect to 
strength, alloy stability, and environmental resistance. 
Certain elements have been found to be deleterious 
to the properties of NiAl, among them are nitrogen, oxy- 
gen, silicon, phosphorus, and sulfur, [jjj The presence of 
these elements must be controlled during the melting pro- 
cesses. On the other hand, some other elements are de- 
sirable because they improve the cohesive properties of 
NiAl. Chromium impurities are important for improv- 
ing its hot corrosion resistance, while boron, carbon, and 
zirconium provide improved resistance of grain bound- 
aries to fracture at elevated temperatures. Jl|,^| Stoloff 
has given an extensive list of atomic additives and their 
effects on the properties of nickel-base superalloys. Q| 

Boron is the main grain-boundary strcngthener in 
NiAl. The beneficial effect of boron additives was first 
observed by Aoki and Izumi in 1979. Boron has a 
strong tendency to segregate to grain boundaries and it 
can improve the tensile ductility of a polycrystal by an 
order of magnitude. This increase in tensile ductility 



is accompanied by a change in the fracture mode from 
brittle intergranular to ductile transgranular |J which 
clearly shows the effectiveness of boron in improving the 
intergranular cohesion in a polycrystal. The strengthen- 
ing effect of boron additives has also been observed when 
they are present as impurity atoms in bulk. The improve- 
ment in strength obtains even with a small concentration 
of boron dopants: 30 weight ppm of boron can give rise 
to a 30% increase in yield strength. || In addition to 
boron, carbon, which is the element next to boron in the 
periodic table, is also a potent strengthener in NiAl. || 
In contrast to boron and carbon, oxygen and nitrogen 
are known to be harmful to the cohesion in NiAl. Indeed 
oxidation is among the most common degradation mech- 
anisms in many metals and alloys. |],[l(| In NiAl, oxygen 
will selectively attack the least noble constituent, which is 
aluminum, and form the stable oxide product AI2O3. pd| ] 
The rate of formation of NiO is negligible compared to 
that of AI2O3. [H§jn|l This strongly-preferential bond- 
ing has also been shown to occur in some recent first- 
principles calculations Jl2|,[l3| and it may be among the 
key microscopic ingredients for the formation of various 
mesoscopic structures (e.g., pores, cracks, and blisters) 
created during an oxygen attack on an intermetallic alloy. 
H In the extreme, oxygen can cause the pesting degra- 
dation phenomenon which happens when some polycrys- 
talline samples are heated in air within a certain range 
of intermediate to high temperatures. Jjjjlyj This pro- 
cess, which is essentially a spontaneous disintegration of 
the polycrystalline alloy to powder, can take place in a 
matter of several hours. 14 15|] 



There have been several first-principles calculations in 
the literature on the effects of impurities on the cohesion 
in nickel aluminides and related alloys. Sun et al. have 
studied the effects of boron and hydrogen on N13AI us- 
ing a full-potential linear muffin-tin orbital (FPLMTO) 
method. Ilq] They emphasized the increase of the inter- 
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FIG. 1. Onsite parameters for orbitals in 16-atom super- 
cells of NiAl with one impurity atom substituting for Al. The 
leftmost column gives the onsite parameters for pure NiAl. 
The lines with dots are the onsite parameters for the 2s, 2p, 
and 3d orbitals at the impurity atom. The horizontal line at 
E F = 1.0475 Ry is the Fermi level of the pure NiAl. The 
Fermi energy for the supercell is 1.0105, 1.0085, 0.9943, and 
0.9927 Ry for B, C, N, and O, respectively. 
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FIG. 2. Site-projected 2s density of states at the impurity 
atom for the case where it is substituting an Al atom. Notice 
the strong resonance at the bottom of the band in the case of 
carbon. The delta functions shown in the spectra for N and O 
are actually very narrow bands (in our supercell calculation) 
with bandwidth of 0.0226 Ry and 0.00198 Ry, respectively. 
The numbers next to the arrows are the partial weights of the 
impurity 2s state in the band, while Et denotes the center of 
the band. The Fermi level of NiAl is at E F = 1.0475 Ry. Each 
spectrum has been given a separate vertical shift for clarity. 



stitial bonding charge as the origin of the beneficial effect 
of boron. Wu et al. calculated the effects of boron and 
phosphorus on the grain-boundary cohesion of iron using 
a full-potential linear augmented plane- wave (FPLAPW) 
method. JlTj] They showed that a combination method of 
the thermodynamic theory of Rice and Wang jl8],[L9) and 
first-principles total-energy calculations can be used to 
determine the grain-boundary embrittlement potency of 
a given impurity. Using the same combination method 
they have also studied the effects of hydrogen and carbon 
impurities in iron and hydrogen, boron, and phosporus 
in nickel. pfj|-p2]| Vacancies and antistructure defects in 
transition-metal aluminides have been studied by several 
different groups. |p^-^5|| 

Previous first-principles studies on the effects of im- 
purity atoms in nickel aluminides have generally focused 
on, and drawn their conclusions from, the calculated total 
energy and electronic charge densities. Insights into the 
bonding and hybridization in the system, however, can 
usually be obtained more clearly by working with local- 
ized basis functions and using the simpler tight-binding 
representation. [p6|-p8| Recently we have shown that ac- 
curate tight-binding parameters can be obtained directly 
from the FPLMTO method. ||§ In this paper we have 
used this method to perform a systematic study of im- 
purities on NiAl. The motivation for carrying out a sys- 
tematic study is the widely different effects that can be 
caused by "nearby" atoms in the periodic table. It is not 
obvious, e.g., why, along the 2p row, boron and carbon 
are good cohesion enhancers in NiAl while the next ele- 
ments, nitrogen and oxygen, are embrittlers. The present 
study has been carried out in an effort to find the answer 
to this question. In the next section, we will give a brief 
description of the FPLMTO method that we use. The 
rest of the paper presents the results of our calculations. 



II. FPLMTO METHOD 

We use the Wills-Price all-electron full-potential imple- 
mentation of the LMTO method. |^|| In FPLMTO, 
no assumption is made about the form of the wave func- 
tions, charge density, or potential. The muffin-tin poten- 
tial is used only to construct the LMTO basis functions 
but the final wave functions, and other quantities derived 
from them, are not limited to such form, p^j Relativistic 
Dirac equations are used for the core states, while the va- 
lence states are treated semirelativistically without spin- 
orbit coupling. For the exchange-correlation potential, 
we use the parametrization of Vosko, Wilk, and Nusair. 
p3[ Within the muffin-tin spheres, lattice harmonics with 
angular momentum I < 8 are used. NiAl is a good para- 
magnetic metal (it has no measurable magnetic ordering 
down to temperatures of a few Kelvin [ j34||35| ] ) therefore 
we do not use spin polarization in our calculation. 
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FIG. 3. Onsite parameters for orbitals in 16-atom super- 
cells of NiAl with one impurity atom substituting for Ni. The 
leftmost column gives the onsite parameters for pure NiAl. 
The lines with dots are the onsite parameters for the 2s, 2p, 
and 3d orbitals at the impurity atom. The horizontal line at 
E F = 1.0475 Ry is the Fermi level of the pure NiAl. The 
Fermi energy for the supercell is 1.0283, 1.0185, 1.0051, and 
1.0014 Ry for B, C, N, and O, respectively. 
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FIG. 4. Site-projected 2s density of states at the impurity 
atom for the case where it is substituting a Ni atom. Notice 
the strong resonance at the bottom of the band in the case of 
carbon. The delta functions shown in the spectra for N and O 
are actually very narrow bands (in our supercell calculation) 
with bandwidth of 0.0225 Ry and 0.00158 Ry, respectively. 
The numbers next to the arrows are the partial weights of the 
impurity 2s state in the band, while Et denotes the center of 
the band. The Fermi level of NiAl is at E F = 1.0475 Ry. Each 
spectrum has been given a separate vertical shift for clarity. 



Impurity is incorporated in our FPLMTO calculations 
by using a 16-atom supercell. |H| NiAl crystallizes in B2 
structure which is a bcc-based structure with one atom 
(Ni or Al) occupying the center of the cube (5,5,5) and 
the other (Al or Ni) at the corner of the cube (0,0,0). 
The cubic supercell is constructed from 2 3 NiAl unit cells 
and the impurity atom is placed at the center of the su- 
percell. Each atom is assigned a minimal basis set con- 
sisting of 9 (spd) orbitals. Since we want to work with 
localized orbitals, the interstitial parameter for each or- 
bital has been uniformly set at k = —0.2 a.u. This gives 
well-localized FPLMTO basis functions with an envelope 
that decays roughly as exp(— \k\v). Q 

The standard FPLMTO method self-consistently cal- 
culates the basis functions, along with the corresponding 
charge density and the resulting total energy, by working 
in momentum space. The program computes the matrix 
elements of the hamiltonian, H a p(k), and the overlap, 
SapCk), matrices from which the energy bands e(k) are 
obtained by diagonalization. From these matrix elements 
in momentum space, we have calculated the matrix ele- 
ments in real space by direct Fourier transform: 



3 



exp (ik • Hj)H a p(R.j). 



(1) 



The onsite parameters are simply the hamiltonian ma- 
trix elements, in real space, between identical orbitals, 
H aa (R = 0). This is computed by using an additional 
code built on top of our FPLMTO program. The distri- 
bution of onsite parameters is an important ingredient in, 
e.g., Anderson's theory of diagonal localization. ]3r^-|39|l 
In this theory, the distribution of onsite parameters, char- 
acterized by the width of the distribution W, competes 
with the strength of the hybridization between the or- 
bitals, which in the impurity case can be taken to be the 
bandwidth of the parent system B. Electron localization 
is more favorable for large values of W/B. p7| 

To obtain a measure of the hybridization strength 
between the orbitals in the system, without having to 
deal explicitly with the multiplicity of hopping and over- 
lap parameters, one can instead examine the density of 
states (DOS) and its atomic-site and angular-momentum 
projections. p5]-^2]| Spin projection is unnecessary since 
NiAl is paramagnetic and we do not use spin polariza- 
tion in our calculations. In this paper the total DOS and 
its projections have been computed using the standard 
tetrahedron method with 35 points in the irreducible 
wedge of the cubic Brillouin zone. The total DOS is 
calculated by summing the contributions from all bands 
and all tetrahedra: [f±3| 



n,k c 



(2) 



where n is the band index while k c is the index for the 
tetrahedra. The site (index i) and angular-momentum 
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FIG. 5. Local charge 
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(in units of the electron charge |e|) 
induced on each atom in the supercell for the case where an 
impurity atom (X = B, C, N, and O) is substitutionary placed 
at an Al site (top panel) and a Ni site (bottom panel). Alu- 
minum charges are marked by filled circles (•), Ni by filled tri- 
angles, while the impurity charges are shown as filled squares. 
In pure NiAl, Al has a charge of +0.2 while, from charge neu- 
trality, Ni has the opposite charge of —0.2. 



(index I) projection of the DOS is obtained by multi- 
plying each contribution with its decomposition weight 
Wmz(k c ) which is obtained from the wavefunctions: 



pu(E) = 2J w ni i(k c ) ■ g n (k c ;E). 



(3) 



Standard DOS decomposition in the FPLMTO method 
differentiates between the muffin-tin (MT) and the inter- 
stitial components of the electron distribution. Fur- 
ther site and angular-momentum decomposition, i.e. the 
calculation of the weights w n i;(k c ), is then performed 
only on the part of the LMTO wavefunction inside the 
MT spheres. The interstitial part is not considered to 
belong to any particular site and therefore is not sub- 
jected to further decomposition. It should be noted that 
this differentiation between MT and interstitial charge is 
an artificial one since it depends on the size of the MT 
sphere which, in common practice, is set rather arbitrar- 
ily by the user of the FPLMTO code. Furthermore, the 
interpretation of such a decomposition is difficult since, 
e.g., the integrated spectral weight for a particular atom 
(vi) is, in general, less than the total number of valence 
electrons assigned to it (rij): 
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W P Z MTO \E)dE< 
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(4) 



Note that the summation over the angular momenta ex- 
tends up to l m , which is a free parameter in an FPLMTO 
calculation (this parameter is set to 8 in this work). In 
general, this parameter is different from (usually much 
greater than) the highest angular momentum L m that 
one uses in defining the FPLMTO basis functions (L m — 
2 for spd basis that we use here). Inside each MT, tails 
from the basis functions centered at other MTs give rise 
to higher angular-momentum harmonics when expanded 
relative to the center of the MT sphere. The parameter 
l m is the cutoff value used in this expansion. p9| ] 

Instead of using this MT decomposition, in this paper 
we have chosen to use an orthogonal decomposition which 
is the one used in tight-binding systems. The FPLMTO 
non-orthogonal matrix elements, i/ Q| a(k) and S a /3(\s), are 
first transformed into an orthogonal system by Lowdin 
transformation. [E3 Since this is a symmetry transforma- 
tion which does not mix components of different angular 
momenta, [ p5p^ ] the weights for the Z-projected DOS 
can be obtained readily from the resulting Lowdin eigen- 
vectors. Details on this scheme have been presented in 
an earlier paper. 0] In this decomposition, the angular 
momentum expansion extends only to L m and the total 
atomic weight is equal to the number of the assigned va- 
lence electrons since the interstitial continuation of each 
FPLMTO basis function has been incorporated properly: 



,(TB) _ 



£/ P T\E)dE = ni 
i=o J - QO 



(5) 



This decomposition method is more appropriate to use 
in our case since we exclusively use localized FPLMTO 
basis functions (specified by negative k parameter). It 
should be pointed out that the Lowdin transformation 
to orthogonal system is used solely to obtain the decom- 
position weight w n u(k c ) for the local DOS; elsewhere in 
this paper we work directly with non-orthogonal TB sys- 
tems. The onsite parameters displayed in Figs. ^ and 
|3|, e.g., are matrix elements of the hamiltonian operator 
in the original non-orthogonal FPLMTO basis; they are 
not, and should not be confused with, the matrix ele- 
ments in a Lowdin orthogonal basis which are nowhere 
presented or analyzed in this paper. 

From the projected DOS, the total number of electrons 
residing on each atomic site can be obtained by integrat- 
ing the corresponding DOS up to the Fermi energy: 



= Pa (E)dE. 



(6) 



Here Ep is the self-consistent Fermi energy calculated for 
each supercell (with impurity atom) and not the Fermi 
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FIG. 6. Onsite parameters for orbitals in 16-atom super- 
cells of NiAl with one impurity atom (Si, P, or S) substituting 
for Al (on the left side of the vertical midline); and with the 
impurity substituting for Ni (right side) . The leftmost column 
gives the onsite parameters for pure NiAl. The lines with dots 
are the onsite parameters for the 3s, 3p, and 3d orbitals at the 
impurity atom. The horizontal line at Ef = 1.0475 Ry is the 
Fermi level of the pure NiAl. In the case where the impurity 
is occupying an Al site, the Fermi energy for the supercell is 
1.0440, 1.0465, and 1.0426 for Si, P, and S, respectively. For 
the case where it is occupying a Ni site, the Fermi energies 
are 1.0765 Ry (Si), 1.0795 Ry (P), and 1.0750 Ry (S). 



energy of the pure NiAl system. In the next section, 
we present the results of our calculations for the onsite 
parameters, projected DOS, and the local charges. 



III. COMPUTATIONAL RESULTS 

A. Boron, Carbon, Nitrogen, and Oxygen 

Fig. [l] displays the calculated onsite parameters in pure 
NiAl and in 16-atom supercells of NiAl with one impu- 
rity atom substituting for Al. We have used the com- 
puted equilibrium lattice constant for NiAl (5.3451 a.u.) 
which is within 2% of the experimental value (5.4450 
a.u.). The XAlyNig supercell (here X stands for the im- 
purity atom) is constructed from 2 3 NiAl unit cells Q 
with the impurity atom placed at (0,0,0); Al atoms at 
(|,0, 0), (i,i,0), {\,\,\), and other equivalent posi- 
tions obtained by permuting the x, y, z coordinates; and 
Ni atoms at (±i,±i,±|) (in units of supercell lattice 
constant). Note that the positions of the Ni atoms are 
all symmetry-equivalent in this supercell. Al atoms, on 



the other hand, occupy three inequivalent sites. This 
gives rise to a small splitting of the Al onsite parameters 
as can be seen in Fig. [|; the corresponding splitting of 
the Ni parameters in the XNiyAlg supercell can be seen 
in Fig. |^. Atomic relaxation has been shown to produce 
only a small change in energy |D| and therefore has been 
ignored in this work. The small value of the computed 
relaxation energy, |l3| and the small size of the onsite- 
parameter splittings in Figs. [I] and provide the justifi- 
cation for our neglect of atomic relaxation in the present 
work. It is unlikely that relaxation will make large quan- 
titative change in, or rearrange the qualitative structure 
of, the onsite-parameter maps in Figs. |l| and [| on which 
we will base much of our discussion in this paper. 

The utility of plotting the onsite parameters system- 
atically, as in Fig. El comes from the fact that it shows 
clearly how well the 2s and 2p parameters of boron match 
those of the corresponding 3s and 3p orbitals of alu- 
minum, and how rapidly this compatibility deteriorates 
as we go from boron to oxygen. To our knowledge this 
almost-perfect compatibility has never been pointed out 
previously in the literature. The onsite parameters for 
the 3d states of the impurity atoms are all much higher 
than the Al-3d parameters. Although results from the 
local density approximation (LDA) for the excited states 
are known in general to be less accurate than the corre- 
sponding results for the occupied states, we believe this 
visible difference is an important feature in explaining the 
efficacy of boron as a cohesion enhancer in NiAl. The 
much higher B-3d parameters would allow the delocal- 
ized B-2s and B-2p states to create wider bands centered 
at their corresponding onsite parameters which, as we 
pointed out previously, match closely to those of Al-3s 
and Al-3p. The overall cohesion is therefore improved by 
increasing the bond order (the difference in occupancy 
between bonding and antibonding states). Gffl 

The bottom of the pure-NiAl bands lies just above the 
zero energy in Fig. |lj. It can therefore be seen clearly 
that the C-2s onsite parameter sits just above this bot- 
tom while those of N-2s and 0-2s orbitals lie below the 
main manifold of pure NiAl. As in standard scattering 
theory, E] this situation opens the possibility for the ex- 
istence of resonance or bound states. In Fig. || we show 
the 2s projected DOS at the impurity atoms. The total 
weight under each curve is equal to 2 (due to spin sum) 
to within 2% accuracy. For this case, where the impu- 
rity atom is occupying an Al site, the nearest neighbors 
of the impurity atom are Ni atoms. The main feature 
of the DOS for B-2s is a broad band which is cleaved 
by its interaction with the neighboring 3d orbitals of Ni. 

This is markedly different from the DOS for C-2s 
in which the dominating feature is the strong resonant 
peak at the bottom of the spectrum. As we move on to 
N-2s and 0-2s, the onsite parameters for these orbitals 
are deep enough to localize the electrons in a bound state. 
This results in a transfer of the spectral weight from the 
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FIG. 7. Site-projected 3s density of states at the impurity 
site (Si, P, or S). A bound state is formed for the cases of 
phosporus atom at Ni site and sulfur atom at Al or Ni site. 
The delta functions shown in the spectra are actually very 
narrow bands (in our supercell calculation) with a bandwidth 
of 47 mRy (P at Ni), 26 mRy (S at Al), and 16 mRy (S at 
Ni). The numbers next to the arrows are the partial weights 
of the impurity 3s state in the corresponding narrow band, 
while Eb denotes the center of the band. The Fermi level of 
NiAl is at Ef = 1.0475 Ry. Each spectrum has been given a 
separate vertical shift for clarity. 



continuum to the bound state. In our supercell calcula- 
tion, the bound state is not manifested as a true delta 
function but it rather appears as a very narrow band 
(with bandwidth of 23 mRy and 2 mRy for N-2s and O- 
2s, respectively, for the case where the impurity is placed 
at an Al site) which is separated by a gap from the main 
spectrum and is displayed as a vertical arrow in Fig. ^. 
This narrow band still contains small hybridization com- 
ponents from other orbitals (this, of course, is just an 
artifact of a supercell calculation) which, as expected, 
diminish as we go from N to O. The total weight of the 
impurity-2s state in the narrow band is displayed next 
to its arrow in Fig. ^ while the rest of the weight still 
remains spread out thinly in the continuum. 

A very similar map of onsite parameters is obtained 
in the alternative case where the impurity atom is sub- 
stituting for a Ni atom, as shown in Fig. ^|. The main 
difference from Fig. [I] is the fact that the impurity-3ri 
levels are pushed to much higher values in this case (by 
about 0.7 — 0.8 Ry). This feature is mainly due to the 
smaller size of Ni, compared to Al, which increases the 



kinetic energy of the orbital (we use a MT radius of 2.30 
a.u. for Al and 1.85 a.u. for Ni). Combined with the 
bond-order argument described previously, this also pro- 
vides a heuristic explanation on why an oxygen impu- 
rity would prefer to occupy a nickel site over an alu- 
minum site, a result which was recently obtained from 
a full-fledged FPLMTO calculation. Q The calculated 
impurity-2s DOS for the case of impurity at Ni site is 
shown in Fig. ^. The general progression from B to O is 
the same as in Fig. ||: A broad band for B-2s, strong reso- 
nance for C-2s, and bound state with an increasing bind- 
ing energy for N-2s and 0-2s. Since the impurity atom is 
surrounded by Al nearest neighbors in this case, instead 
of Ni atoms, we do not see as strong a band cleavage 
around the energy of the Ni-3e? orbitals as seen in Fig. g. 
Below the Fermi energy, the similarity of the results ob- 
tained for impurity at Al and Ni sites shows that these 
features, e.g. the compatibility of the onsite parameters 
for B and Al in NiAl environment, are largely indepen- 
dent of the atomic arrangement in the crystal. This is 
not surprising since an onsite parameter is sensitive only 
to the average potential at its atomic site. This suggests 
that our results in this paper, which have been obtained 
for NiAl host using 16-atom supercells, may have some 
relevance also to other nickel-aluminide alloys with dif- 
ferent concentrations of impurity atoms. 

Fig. H shows the charges induced on each atom in the 
supercell which have been calculated by substracting the 
total number of electrons on the site, qi in Eq.(^|), from 
the assigned number of valence electrons rij : 



Qi- 



(7) 



In all cases, Ni is found to be more electronegative than 
Al (the Pauling electronegativity of Ni and Al is 1.91 and 
1.61 respectively [^0j). In pure NiAl, Al has a charge of 
+0.2 (in units of electron charge |e|) while, from charge 
neutrality, Ni has the opposite charge of —0.2. The im- 
purities from the 2p row that we have studied in this 
work have Pauling electronegativity of 2.04, 2.55, 3.04, 
and 3.44 for B, C, N, and O, respectively. |5(J It can be 
seen that this electronegativity trend is followed rather 
well in Fig. |[ In the case of impurity atom at Al site 
(top panel in Fig. ||) , N and O are sufficiently electroneg- 
ative to change the sign of the induced charge on their Ni 
nearest neighbors, relative to the sign of the correspond- 
ing charge when B or C is present. Thus a portion of 
the valence electrons localized at the N or O bound state 
comes from their nearest-neighbor Ni atoms. In the alter- 
native case where the impurity is occupying the Ni site, 
a jump in the induced charge on the Ni atoms is clearly 
seen in the bottom panel of Fig. ||. Although in this case 
they are no longer the nearest neighbors of the impu- 
rity (since they are separated from it by the Al atoms), 
the formation of the bound state in N and O still has a 
substantial effect on the Ni atoms. Two reasons may be 
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given to explain this strong interaction between the im- 
purity and the Ni atoms. First, the Al nearest neighbors 
are already positively charged, therefore it is relatively 
harder for the impurity atom to attract their electrons. 
Second, the DOS of NiAl is dominated by strong Ni-3d 
peaks which are situated just below the Fermi energy. 
p3[ These peaks are sufficiently wide to suggest that the 
Ni-3d electrons in this alloy are well delocalized. Their 
proximity to the Fermi energy then strongly expose them 
to changes in the potential as that caused by the forma- 
tion of a bound state on a nearby atom. 

B. Silicon, Phosphorus, and Sulfur 

The elements from the 3p row of the periodic table: 
Si, P, and S, have been known to be strong embrittlers 
in NiAl. jjj It is therefore interesting to examine whether 
the correlation that we have obtained in the previous sub- 
section between the matching of the onsite parameters 
and the macroscopic embrittling/strengthening potency 
of the impurity persists also for these elements. Fig. 
displays the calculated onsite parameters for orbitals in 
16-atom supercells of NiAl containing one impurity atom 
(Si, P, or S) which substitutes for an Al (shown on the 
left side of the vertical midline in Fig. 0) or a Ni atom 
(shown on the right side of the midline). A major differ- 
ence from the corresponding plots of onsite parameters in 
Fig. and Fig. is the decreasing trend of the impurity- 
3d levels as we go to higher atomic number (from Si to 
S). This is due to the fact that the basis orbitals that 
we use in this case (3s, 3p, and 3d) all have the same 
principal quantum number. The difference in their lev- 
els therefore originates mainly from the difference in the 
effective centrifugal potential (the 1(1 + l)r~ 2 term in 
the radial Schrodinger equation) , which is independent of 
the atomic number. |5i| , |5^] In contrast, the basis orbitals 
that we use for the 2p elements in Fig. and Fig. (2s, 
2p, and 3d) come from two different principal quantum 
number shells. In this case, in addition to the centrifugal 
potential, the splitting among the onsite levels is also de- 
termined by the Coulomb potential of the nucleus which 
increases with the atomic number. Thus the 2p level de- 
creases in concert with the 2s level while the splitting 
between them and the 3c? level increases with the atomic 
number as we go from boron to oxygen in Figs. and |. 

Fig. shows the resulting projected DOS for the 
lowest-lying valence (3s) state of Si, P, and S at the im- 
purity site. As in Fig. |2| when the impurity is placed at 
the Al site, its DOS features a peak close to the Fermi 
energy due to its strong hybridization with the 3d states 
of its neighboring Ni atoms. Except for the case of Si at 
Al site, where the resonance at the bottom of the spec- 
trum is relatively weak, the DOS curves in Fig. are all 
dominated either by a very strong resonance (P at Al; Si 
at Ni) or a bound state that is completely separated from 



the main spectrum (S at Al; P at Ni; and S at Ni). Sil- 
icon, phosporus, and sulfur are known to be embrittlers 
in NiAl. JjJ These results therefore support the correla- 
tion that we have obtained in the previous section that 
relates the localization of the valence electrons at the im- 
purity site with the macroscopic embrittling character of 
the impurity atom when it is present in NiAl. 

Although the weak resonance in the case of Si at Al 
site seems to defy this correlation (note that, for reason 
of presentation clarity, the projected-DOS curves that 
we show in Figs. ||, 0, and have been obtained by 
convoluting the FPLMTO DOS with a Gaussian smear- 
ing function of width about 10 mRy), it should also be 
noted that its 3d-state level in Fig. |6| is much lower than 
the corresponding 3d level for, e.g., boron or carbon in 
Fig. 0. As has been pointed out in the previous sub- 
section, this much-lower 3d level exerts an 'onsite pres- 
sure' on its lower-lying s and p states against forming a 
wider band (due to its orthogonality with these states). 
This results in narrower bands under the Fermi level and, 
consequently, in reduced bond order and weaker metallic 
character of bonding around the impurity site. This may 
explain why silicon is an embrittler in NiAl while car- 
bon, which has a similar set of onsite parameters below 
the Fermi energy as can be seen by comparing Fig. and 
Fig. 0, is in contrast a cohesion enhancer. 

IV. SUMMARY 

In this paper we have performed a systematic study of 
impurity hybridization in the refractory alloy NiAl. Im- 
purity atoms from the 2p row (B, C, N, and O) and the 
3p row (Si, P, and S) of the periodic table have been ex- 
amined. The purpose of this study is to understand the 
origin of the embrittling/strengthening property of impu- 
rity atoms in alloys in terms of the compatibility of their 
onsite parameters and their orbital hybridization. We 
found that the onsite parameters of boron, which is the 
prime cohesion enhancer in NiAl, are highly compatible 
with those of the NiAl host below the Fermi energy. In 
addition, its higher-lying atomic levels are located higher 
than the corresponding levels for Al. This allows the 2s 
and 2p states of boron to hybridize more strongly with 
the orbitals at the neighboring atoms, form wider valence 
bands centered below the Fermi energy, and increase the 
bond order. These two properties, the compatibility of 
the onsite parameters and the relative location of the 
higher- lying states of the impurity atom, have been found 
useful in understanding the electronic structure of the 
impurities and their effects on the cohesion in NiAl. 
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